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A ♦. heoi'ii  t i ca  I ar.aly<*in  is  present  id 
for  the  jr'as  i f icat  ion  process  of  condensed 
nateria?  suddenly  exposed  to  a hot  stap- 
n a n t pas. 
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When  the  pas  i fica* Ion  rate  is  piven 
by  an  Ai'rhenius  j -.-w  the  surface  terpera- 
ture  rises  to  th^  non - pa ? : f i c a t i on  jump 
value  upon  contact  with  thf  hot  pesos.  The 
convective  effect  due  to  surface  ^ep^c•s- 
s i o n and  p a s i f i c a l i o n flow,  plus  the  o n e r - 
gy  absoi'b»^d  bv  the  pas  ; f ica  t ier.  causes  the 
surface  t empera t ure  to  dccrc  i* e with  time. 
An  analytical  description  of  thv  surface 
temperature  aiK*’  cor  cent  I’at  i on  history  is 
presented,  topp'th.er  with  closed  form  ex- 
pressions for  thi-  temperature  and  concen- 
tration profiles  for  small  atai  large  times. 

A closed  form  analytical  description 
is  also  piven  for  the  pasification  proc<^r.s 
tarV-r  surface  c qu  1 1 i t r i uia  "Clasius  Clapev- 
ron”  cenJitioriS.  In  this  cane  the  5-urface 
temperjture  ’■•c::ains  at  a constant  jump  va^ 
i.e,  and  the  tetiipoi  at  ure  and  concentration 
profiles  h2v»'  a rinllarity  form. 
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I , Introduct ion 


During  the  past  ycar*s  a considerable 
effort  has  been  devoted  towards  uii  J rs  t jtM 
ing  the  ignition  mechanism  of  a condensed 
material.  A large  amount  of  data  has  been 
gathered  e xpe r imen f a 1 1 y and  seveial  theo- 
ries have  been  propoiacd  to  ir.tcr-prt't  these 
data.  Un f ort una t e 1 y , however,  ignition  is 
a very  complex  phenom^*non  and  many  d : f f r*. r*t 
physical  and  cheniical  processes  dtx  involv* 
during  an  ignition  c xpor  i me n t , The  analyt- 
ical tools  available  are  not  suited  to 
tackle  the  complete  prol'lcrr.  aii  i several 
simplifications  have  to  be*  introduced.  In 
addition,  some  of  the  physical  and  chemi- 
cal paranoters  involved  in  the  ignition 
event  are  not  know*  i li  a given  experimen- 
tal setup  atjJ  this  picvents  meaningful 
comparisons  witfi  r- x ^ .;.r i ti.o r t t 1 data.  Actual 
ly,  thest'  data  arc-  cl  ten  used  to  derive 
the  unknown  parameters,  thereby  obtaining 
good  ag.r^'*r?nt  with  those  particular  data. 

At-  the  present  ♦^ime,  no  consonstjs 
exists  as  to  the  location  of  the  reaction 
respr.nsiMe  for  ignition.  In  fact,  one  of 
the  main  oljectivcs  of  both  t h oo re t i c i a ns 
and  exper i men t a 1 i s I s working  on  ignition 
reseai'ch  has  been  to  elucidate  where  foi 
a class  of  propellants  igniti.'H  tales 
place.  To  achieve  this  goal  it  is  neces- 
sary to  develop  an.ilysi.s  with  the  same  a:- 
sumptions  but  differing  only  on  the  loca- 
tion of  the  key  reaction  leading  to  igni- 
tion. Tt  appears  that  a p j'l-ox  i 1 1 c anal- 
yses yield  a better  un  lerst  and  i tsg  of  th*' 
phenovn-notj  under  study  thin  a l.u'go  amount 
of  numcri-'al  r»rultr.  , although  the  accui 
cy  of  the  lattCT'  for  a given  s».  t of  pir.*- 
mefers  is  superior  to  !bal  of  tlie  approxi- 
mate methods.  Howrvri-,  asymptotic  tech- 
niques hav(  been  (!<vcloped  to  a point  su*  h 
that  while  posuessing  good  aec  ut  tev  ( h<  v 
still  can  retain  as  g.ener'Jl  as.sur.pt  ic'n?  r 
those  useil  in  the  usu.O  numeric  an-il* 
yscs.  These  a s n u-.p  t i on  s shvnjld  b«.  as  j-e,il 
istic  nr  possible  .ind  yrt  simple  enough  i '■ 
be  trc.ital'le  bv  a.ympioiic  t e.  hj>  iques  , ."'n  1 
f u r i.  h<’ rmo  r«»  they  1 »*  .imeniMt-  t c*  I'C 

e xi'C r in.e n t a 1 1 v r-e pr  od uc e d . 
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ihc-  re  s a 1 1 .'J  ' o f th^*3r  a •^yir.j-tot  ic  anal- 
yi.oii  should  b'.*  chocked  vith  r.ood  r.ur»**  r i v.  a \ 
solutions  ai.d  then  ar.  - xpe  r i no t , r • p r’Odu  ^ 
as  closely  as  possIMe  th*'*  ansurpt  iens  of 
the  analyst’s  should  hv  carried  oat  to  ob- 
tain not  only  th*.-  if^nition  tine  but  is  rr.ariy 
sifnificant  viriallos  (such  as  r o npo  r 1 1 u r-e 
and  r.ass  fructio.is  profiles')  as  possible. 
Only  then  jne  an  i np,  f u I experimental  ron^cir: 
sons  will  be  possible  in  trying,  to  aset-r- 
tain  the  location  of  the  reaction  producim’ 
i g n i t ion. 

The  ignition  even*  starts  with  the  ap- 
plication of  li  ntin’.jliic  that  rnhjnces  the 
gasification  of  tb  material  anl  a s<-ltac- 
c e 1 e ra  t or y e xo  t h**r  n i c che n i c a 1 p r'oee  s s . 
Asy.if  t o t i c analyj-s  h.v»?  be^n  carried  out 
assuming  t h ;?  sane  ignition  stimulus,  name- 
ly a constant  flux  of  radintil  etjorgy  ab- 
sorbed it  the  sur^'ace  but  co  nr.  i de  r i n g that 
either  a condensed  phase*,  heterogeneous^ 
or  gas -phase  ^ ^ ^ exothc'-mic  re.iCtion  is 
rcspotisible  for  ignition.  We  are  develop- 
ing similar  analyses  for  an  ignition  stim- 
ulus consisting  in  the  iastar. t<. neous  con- 
tact of  the  cotidc’nsed  material  with  ct  hot 
gas.  Tais  ignition  stimulus  can  be  rep'ro- 
ducod  e xper i men t a 1 1 y by  reflecting  a shock 
wave  on  a condensed  material  located  tat  the 
end  of  a shock  tube.  Interest  in  studying 
this  j)articular  type  of  ignition  stems  f ron 
the  c on s i df rab 1 c amount  of  existing  experi 
men  t s ^ ‘ ^ , nuinerical  results' arid  ap- 
pro x i m ■»  f.  e e.naly ses''^  » ' ' * ' ^ ^ that  will  fa- 

cilitate the  task  of  dcr  ivinii,  couclusious 
regai’diug  the  ignition  mechanism  in  a par- 
ticular material.  We  have  carried  out  re- 
cently asyrrptotic  analyses  for  large  acti- 
vation energies  of  the  ?;0 1 ei'oge  neou  r>  ^ and 
condensed  phare^'  shock- tube  i gn  i t ion  , v;  i t h 
the.  idea  of  generating  simple  analytical 
descriptions  of  the  ignition  process,  and 
the  present  p-aper  is  a first  step)  in  a 
p'arallel  study  of  gas-nhase  ignit.ion.  As 
pointed  out  previous! y^  , gas-phase  ignition 
Ificory  is  mere  co:n})l  ic.ated  than  either  the 
heterogeneous  or  condensed  phase  theories, 
since  it  is  necessary  to  consider  tv;o  dif- 
ferent chemical  processes;  gasification  and 
ignition.  The  ignition  event  vrill  be  con- 
trolled by  either'  of  these  sequential  pro£ 
esecs  depending  on  their  relative  rates. 

The  analysis  of  the  gasification  process 
reveals  several  intererting  features  which 
mak'-s  it  advantageous  to  present  separate- 
ly such  an  analysis,  to  be  followed  by  the 
study  of  the  gas  p^ha.se  ignition  process. 

The  gasificition  reaction  may  be  consid 
tred  as  s rate  process  or  as  an  equilibrium 
pu'cccss.  Usually  'equilibrium  v^ill  be  main 
tained  only  in  the  case  of  liquids.  The 
analysis  in  the  main  text  assumes  a rate 
process  described  by  a 7:ero  order  endother 
flic  A»''rhenins  rate  lavr.  Mov/evei'  in  an  ap- 
p>endix  we  develop  n par-allel  analysis  for 
the  eq»ii  librium  case,  wlien  the  Clasius-Cla 
peyron  equation  is  used  to  r’clatc  the  sur- 
face mass  fraction  of  gasified  liquid  to 
the  surface  temperature . 

In  the  analysis  we  take  inio  account 


the  conve-'tive  cfi^'cts  clue  to  the  flow  jr. 
suiting  from  the  gasification  ‘A  the  cori- 
densed  material,  torcther  v/Ith  the  unst"  dy 
effects  of  h»'jt  cof.duction  in  the  solid  and 
rrs  phar. es  a<id  diffjsiori  in  the  gas  phace. 
Tiic  conser  V-’ t ion  equitions  are  .solved  by 
asymptotic  t ec  ht;  i qu*‘ based  on  th<*  re  jH  stir 
assun.pil  ion  that  the  rior.  d i n'' ns  i ona  1 activa- 
t ion-en'‘r'’.y  of  the  gasification  reaction, 
or  r-atio  of  activation  tempei-uturc  to  t).*- 
solid  initial  tci.perat uro , is  large  com- 
pare'll to  unity.  In  the  context  of  this 
asvmp'totic  analysis  a.n  int'^gral  equation  is 
derive!  wliose  solution  yields  the  surface 
t ? n.p^  ra  t i' : e as  a function  of  time;  the  sue 
face  t c‘ n.j'f  !•  j t u re  rises  instantaneously  at 
t = 0 to  a jump  value,  and  then  decreases 
slowly  as  a result  of  the  endothermic  gas! 
fication  process  which  ?n  addition  blows 
the  hot  gases  away  from  the  solid  surface 
thus  reducing  the  heat  conduction  flux. 
Asymptotic  exprcssioris  giving  the  tempera- 
ture and  r..ass  fraction  profiles  for  small 
and  large  values  of  lime  are  also  given  in 
the  paper. 


II  . 


Formnl at  ion 


We  consider'  the  case  in  which  a con- 
densed mjicrial  at  a temperature  T^,  occu- 
pies the  half  space  x . At  tiu.e  r.ero  the 
t empe ra t u re  of  the  addnccnl  stagnant  gas 
that  occupies  t*ne  half  space  x >0  is  sud- 
denly raised  to  a high  value  Tp^,  to  pro- 
duce a cor',  espoud  iup,  rise  in  surface  tempe 
ratui’o,  and  thereby  cntiancc  a g a s i f i a t .i  on 
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rig.l  Schematic  representation  of  the 
gasification  proccvss. 

as  inert.  For  convenience,  the  origin  x -0 
i r,  maintained  at  the  rurface  of  t))0  re- 
gi'cssing  materiYl  and  the  mass  coordin.ilo 
4'  is  introduced  in  canalyzing  the  gas  phase 

We  will  use  the  ’••11  jurtifird  assump- 
tions that  the  pressure  is  cotislant  and 
that  th»’  ''ffect  of  exi<'rr.al  foT'cer  and 
viscous  stresses  arc  ticflfgible.  We  use 
Kick's  } w to  calculate  the  diffusion  ve- 
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xocitiow  and  consider  iho  specific  heats  to 
be  constant  in  eacli  phase.  It  is  also  as- 
suired  that  the  density  and  tliermal  conduc- 
tivity in  the  condensed  phase  as  veil  as 
their  product  in  the  p.as  phase  are  constant. 
Under  these  conditions  the  conservation 
equations  foi'  energy  in  the  solid  and  gas 
phases,  and  species  conservation  equations 
for  oxidizer  and  fuel,  are  respectively 
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An  equation  identical  tc  (3)  exists  for 
the  oxidizer  mass  fraction  , except  that 
the  oxidizei'  diffusivity  should  replace 

Dg  . The  initial  and  boundary  conditions 
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v.'here  = Bexp(-E/RTj)  is  the  value  of  v 
for  , equations  (l)-(8)  become 
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is  ‘ g g <* S’  I j, 

where  all  the  symbols  are  defined  in  the 
nomenclature.  The  sui*face  regression  rate 
V is  given  by 

v = Fiexp(-K/KTg)  ( 9 ) 

and  b is  the  heat  of  gasification  referred 
to  the  "junp’'  tempei-at  ure  Ti,  which  is  the 
one  established  at  the  sui’face  immediately 
after  the  application  of  the  hot  gases.  In 
absence  of  ga:,  i f ica  t ior.  the  surface  tempe- 
rature remains  fixed  at  its  "jump"  value 


7,  = T + It  -t  ] (10) 

1 o '•goo-'i  + r 

In  the  analysis  that  follows,  the  quan 
tity  Pj,Xg  is  taken  an  constant,  and  its 
value  IS  chosen  in  such  a way  that  the 
correct  "jump'  temperature  (with  arbitrary 
dependence  of  on  temperature)  calcu- 

lated in  reference  ?0  is  recovered,  when 
the  inert  problem  is  solved  in  the  assump- 
tion of  constant  thermal  nropertios. 

Introducing  the  nond  imensionjl  varia- 
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where  v/v^  r exp  (-E/RT^  + b/RTj)  . 

An  equation  identical  to  (16)  defines 
the  oxidizer  mass-fraction  , provided 

that  the  Lewis  number  for  the  oxidizer  be 


is  used  instead  of  be. 


III.  Surftice  Tempera tiire  and  Regresrion 
Rate  M -story. 

At  small  times  the  temperature  gratli- 
ents  at  the  surface  are  large  and  there- 
fore the  energy  absorbed  by  the  gasifica- 
tion reaction  is  negligible  compared  to 
the  energy  transport  by  conduction  from 
the  gas  phase  to  the  condensed  plwise.  Upon 
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canlaci  with  the  he:  ['rtJ.'C-s  tho  suiface  tvm 
por^»lui't-  will  be  i p.M an t an ■r'o us  1 y raised  to 
a ju'r.p  VmIuo  correspond  i to  z^-io  p.asifi- 
cation  rate  hu  t , s ubso  q uo  :i  i ] y , r,ar>  i f i ca  t ion 
introduces  two  effects  that  cause  the  sur- 
face t orrpcr'a  t u i‘v  ind  the  re;»rossion  rate 
to  decrease  with  tinie.  Tir-st,  a fraction 
of  the  enerr,y  cominp  iron  the  pas  phase  is 
absorbed  by  the  gasification  reaction,  so 
that  the  surface  temperutui'e  can  not  be 
nioitrtainod  at  its  jump  value.  Second,  the 
convective  flow  resiiltinp  from  the  con- 
densed phase  gasification  blows  ttic  hot 
gas  away  fror  the  surface  and  brinp.s  cold 
mat»>rial  towai’ds  the  surface,  thei>eby  re- 
ducinp  the  surface  temperature  and  the  re- 
gression rate. 


1 n t ro  due  i n g the  non  - d i ne  n s i on  >i  1 varii 
bles  (?3)-i?hJ  ar.d  taking,  the  liv.it  f -*■  « , 
equationr  (1m  )-(?!)  simjlify  to: 
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3 c IiC  3 z ^ 
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x(o,0)  = x^ia  ,0) 

( 29  ) 

We  shall  analyrc  first  the  period  of 
time  when  the  surface  regression  rate  is 
r’ed’iced  by  a factor'  of  or-der  unity  from  its 
’’initial"  value  corresponding  to  the  jump 
temperature  Ti . For  large  values  of  the 
non-d  imens i onal  activation  energy  F/KTi  , 
small  increments  in  the  surface  temperature 
of  the  oi-dor  RT^/L  pi'oduno  a change  in  the 
gasification  ra^e  by  a factor  c . For  those 
early  times  the  corresponding  deviations 

of  e and  0 from  their  inert  values 

g 

3^  = -0,/r  = c:'f(f./2vV)  (22) 

are  sr.all,  of  order  r.'l' ^ / { g T j ) } , so 

that  it  is  convenient  to  describe  this  ear 
ly  period  in  terms  of  the  non - d ime n s i ona 1 
V a r i a b 1 c s 

X = 6(0^-0)/r  , Xg 

y = 6Yp/[arOTD/L7]  (24) 


0 = , z = BC/r  (21) 

whei'c 

B = e(Tj-T^)/RTj  (26) 

is  a non-dimensional  ictivation  energy. 

The  now  non  - d i mens  i c:'.  a 1 lime  o and  fuel 
mass  fraction  y have  been  defined  so  that 
they  are  of  order  unity  during  this  period. 
The  analysis  that  foJlov;s  is  based  on  the 
assumption  that  the  par.imetcr  6 is  large 
compared  to  unity  which  is  equivalent  to 
require  I h 1 1 the  gas  i f ica  i.  ion  rate  at  the 
jump  temperature  be  much  faster  th<ifi  the 
gasification  rate  at  the  initial  tempera- 
ture. Thus,  it  is  necessary  to  aesumo  that 
6 is  large  in  order  to  hive  a physically 
realistic  problem  in  v^hich  a non  reacting 
condensed  material  i r.  suddenly  br'ought  to 
vigoi'ous  g 1 r.  i f i ca  t i on  upon  contact  v.ith  hot 
gases.  The  ratio  of  thermal  i*e  s pc  n s i v i t ies 
r is  usually  small  , but  it  may  bcconio  of 
order  unity  for  large  valuer,  of  pressure 
or  T , since  it  is  propoi’ t i o na  1 to  the 
square  root  of  p A . Thus,  the  par.irneter 
p/r  is  .:it  least  of  order  6 . 


+ r 


( J 4 r )c>:[i  ( - Xs  ) 


( 3C  ) 


lx 

3z 


s 


/lTc  exp("Xj.) 


(31) 


The  factor  exp  (-Xr.)  ” v/v-  account'-  for 
the  reduction  of  gasification  rat"  with 
t empera  t ur o . 

The  parameter  T appe-iring  in  i-quation 
(27)  for'  tlv’  case  of  condensed  phase  temp£ 
rature  x «*quals  ar(l-il  ).  Vho  < xintjon  for 
the  gas  phase  ter.perat  utw-  Xp  i idv.-tical 
to  Eq,(?7)  with  P replaced  by  oifl+I')/!. 


Only  the  gradient  of  the  inert  tempe- 
rature appears  in  the  effi'-*  of  convection 
given  by  the  I'ight-har.d  side  of  Fq  (27), 
and  no  convective  effects  appear  in  bq  (20) 
describing  tlic  mass  fractions.  It  sViould 
also  be  observed  that  in  the  limit  £ -■  ® 
the  equations  and  boundary  conditions  de- 
scribing the  temperature  d i s t r i I'li  t i on  are 
uncoupled  ircin  those  describing  icasx  frac- 
tion. Green  functions  may  be  us*>d  to 
write  the  solution  to  equations  (27),  (29) 
and  (30)  as 


xV  ,z) 


p_ 

2ti' 


('  cxp(-x')  ” 

S f CXIl(-Z  '^/Mo' ) 

0 /a  ' ( 3-o')  p 


X G(  7 ,c  r.  ',o')dz'dc'  + 


a e X p f - z •’  / 4 ( o - o ' )^ 

0 /ti  ( 0 - Cl  ' ) 


( 32  ) 


wl.ore 

G ( 7. , cl ; :c ' , ')  ' ) = e X p f - ( z - z ' ) ’ / 4 ( o • n ' ^ •( 

4 o X p f ( z + z ' 'V  / 4 ( 0 - o ' )J  (33) 


f n c r; 


I’Vci ) 11  1 1 i np,  •’(ju.ition  (32)  .i1  ttiP  r.ut- 
iiid  uiiiiip,  pcjuatioii  (30)  o:>  imcj-.i’''! 
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ir.  IS  equ.ition  aescrjDcs  me  evoiui  ion 
time  of  the  surface  temijrrature  as  a f 
t ion  of  the  unique  pai’aneter  A - a(r  + r 
which  is  a ncasure  of  tiio  coolin^»  effe 
the  pasificaticn  flow.  The  second  ter 
the  right  hand  side  of  equation  ( ) a 
counts  for  th?  heat  dhsorbed  Ly  the  en 
thermic  gasification  and  vfill  bo  the  o 
one  appearing  if  the  convective  effect 
were  neglected. 

From  equations  (2d),  (29)  and  (31) 
integral  equation  may  bo  derived  relat 
the  surface  mass  fraction  to  time,  vi? 

0 exp(-x^) 

yg  ■ / - — : ^ Jc  ' 

o /flto  - 0 ' ) 


This 
late  the 
fraction 
y in  t h 

r 

observed 
suT-f  ace 
in  the  r 
with  t li  o 
the  1 i rn  :■ 
are  n e p.  .* 


equation  can  al sc  be  used  to 
surface  value  of  the  oxidizetr 
, if  we  use  1 - instead 

e definition  of  y (Eq.24).  Jt 
that  the  mass  fractions  at  th 
cire  proportional  to  the  second 
ight  hand  side  of  (34)  associa 
gasification  reaction,  since 
t 6 +•  » , the  convective  c f f e c I 
i p.  i b 1 e , 


A fin 
ns  (34 
egrate 
t used 
h an  i 
and  ( 
“d  imen 


with  the 

rl  S C S t 


difference  versi 
d (35)  fi  a G been 
rougli  <1  prv')cedur 
re  f e r*.  nee  1 . Th 
ration  ’ re  slio wn 
or  several  value 
al  regression  va 
plotted  in  fig 
urn  a I t ho  t i mo  o 
gas  and  cont  f nuo 
d s zero  a f t erwar 


o one  o t t 


jTosent  study  from  its  ii  r p 1 i ca  t ici. 

to  p<*.  r phase  i^,nition»  it  is  i !*  t o res  t i n 
to  ohst. rva  that  under  cei'tiin  pai'ticular 
conditions,  the  present  analysis  enn  be 
easily  ^»enorilized  to  account  for  fast  gas 
phase  reart  Ions  and  even  to  describe  the 
ip.nition  event,  bimely.  If  tlie  fuel  react'’ 
immediately  a f 1 1:  r'  gasification  w i t h t h e 
ambient  oxidi.’.er,  the  heat  rel^.ured  by  this 
exothermic  reaction  shcrld  be  subtracted 
from  the  heat  absorbed  by  t h -•  p.i  s i f ica  t i o n 
reaction  L.  It  is  not  difficult  to  shoe 
that  in  this  case  Yp-0  in  the  gar.  phase 
and  the  ratio  ( v ) , where  v is 

the  mass  stoichiometric  ratio  o x i d i /or/f  u el, 
satisfies  hq  (3)  and  the  sar.ie  initial  and 
boundary  condition  than  the  fuel  mass  fi'ac 
tion  in  the  non-reactivc  case.  In  par-tic- 
ular  if  the  net  value  of  L is  positive 
rio.3  giving  used  to  calcu- 
late Y (a)  if  wo  write  (Y  -Y_  )/(Y  +v)= 

Os  0»'0:5  O” 

= y ft  r ( i + T’ ) /Lo /C  . The  arialysis  ceases  to 

be  valid  if  the  resvlting  value  Y^^  become^ 
negative,  when  the  reaction  zonfj  moves 
away  from  the  surface.  J 

for  sufficiently  cxothvrnic  reactions 
the  net  value  of  h is  negative;  under 
these  conditions,  A is  negative  and  the  sun 
face  temperature  satisfies  Kq  .(34)  if  its 
left  hand  side  is  replaced  by  (-x  This 

equation  has  bo:*!k  numci'ieal  ly  integrated 

in  reference  20, tc  generate  0.  as  a func 

1 g n ““ 

tion  of  A.  It  is  observed  that  ignition 
occurs  only  for  values  cf  -A  < 2,  in  which 
case  “X^.  Is  positive  and  ^-.rows  to  infinite 

at  a finite  time  0.  (A). 

ign 

^ * Z.®  utu  ond  Concentration  Profi  Ic  s 

The  temper' a turn  and  concenti'ation  pro- 
files can  be  calculated,  ir\  principle,  by 
using  Eq.(32)  in  terms  of  the  calculated 
surface  temperature  history.  Hcv;ever,  we 
shall  instead,  derive  expansions  for  tlie 
tempers t ure  and  concentration  profiles  for 
small  and  large  values  of  the  non  - d irnens  i on 
al  time  0.  In  tliese  expansions  the  non-di 
mcnsional  distance  to  the  solid  surface 
is  scaled  with  the  characteristic  thick- 
ness of  the  heat  conduction  layer  2 /o  t o 
pZ'oduce  the  similai’ity  variable  n = 7/2/a" 
for  staall  and  large  times.  shall  later 

use  these  expansions  to  generate  approxi- 
mate expressions  that  can  be  used  to  calc£ 
late  the  surface  t oir.peia  t ur  c for  all  tines. 

For  small  values  of  o the  factor 
cxp(-Xg)  Fq.(27),  accounting  for  the 
change  in  regression  velocity  v’ith  chang.es 
in  surface  temperature,  can  be  written  e- 
qual  to  1 and,  in  addition,  the  boat  of 
gasification  isi  Kq.(30)  can  be  neglected; 
the  resulting  syst*'m  of  ^'quat  ions  has  a sim 
ilarity  solution  of  the  form  x = th)  , 

where 

= [(  2 + A ) //irj  c xp( -n^  ) -t  ( 2- A ) ner  f cn  (37) 


In  order  to  oltain  a more  accui-ate  expres- 
sion cf  X ioi'  small  o v;e  can  use  the  expan 
s ion  ~ 


X = (38) 


which,  when  sulstituted  in  Lqs.{27),  (29) 
and  (33),  previously  v;i*itten  in  terms  of 
the  variables  o and  fi  leads  to  a system  of 
e q a t i o r.  s , v;  i t h solutions  given  by  E c s . ( 3 7 ) 
and 


L-2-A 


4 ( A -*■  2 ) ( F - A ) I For  fen 
^ 3tf  J L 2 


n ^ e r f c n 


. ^ (4.2)oxp(-n2)  (39) 

where  P equals  ar(l-ir)  in  the  condensed 
phase  and  ar(l+r)/r  in  the  gas  phase. 


Thus,  for  early  times  the  surface  lera- 
perature  is  given  by 


Xs  =\/f-  (2.M-a[l  1 f t |A  (A. 2)]  . --(.0) 

This  expansion  can  also  be  derived  dii-cctly 
from  equation  (3h).  It  is  plotted  for  com 
pai'ison  with  the  exact  solution  in  fig- 
ure (2)  for  the  case  A - 0.  Similarly , the 
early  times  mass  fraction  profilefi  are 
given  by 


2.'^ 


- porfc(w) 
/tt 


• 0 ( 2 ■(  A ) r»  * t r J c ( p ) - exp(-vi^)  + 

^ /7 


erfc(  VI ) 


>-]* 


(41) 


whore  p = n /Lc . 

For  large  a the  surface  value  of  x 

is  large  and  it  becomes  a slov'ly  vai-ying 
function  of  o . If  Fqs*  (27)  and  (28)  are 
written  in  terms  of  the  similarity  varia- 
ble n and  0 , the  time  derivative  term  in 

the  equations  can  be  neglected  in  first  up 
proxihuition  for  large  times,  and  the  tempo 
raturc  and  concentra t ion  profiles  take  a 
quas  5 - £ iri5  i 1 a r i t y form.  A more  accurate  rep_ 
rcsentaticvi  of  tlie  temperature  and  conct". - 
tratioii  profiles  can  be  obtained  as  expan- 
sions of  X » X -tid  y in  dccnrisinr  powers 
of  X times  fun’etions  of  the  variable  n , 
Thus  when  the  expa!>5ion 


X " x^/u)F^(n)  + r^(n)  + (i/x^)*'o*h)  + — 

(42) 

and  a similar  expansion  in  powers  of  Xg 
fo!'  Xg  introduced  Jn  Fq.(30),  th^' 
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re]  a t io;: 
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/S' 


is  c t a i n e d , w h c* r- e 


r'.(o)  i IT’,  (o) 

^ ±i- 

2 (!  1 + T5 


■;  '1 3) 
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and  the  subscript  g indicates  gas  phase 
variable.  When  bqs(27)  and  (29)  a!e  also 
written  in  terms  of  r,  and  Xg  ‘■''iid  x ‘’■iid  Xj. 
are  sul  s t i t ut«l  by  their  expansion  ('*2)  we 
obtain,  by  collecting  lil>e  powers  of  '7 

system  of  equatioiiS  with  the  sclutions 

= crfc(n)  + a^  [^e  xp(  - ) - cr  f c ( >1  )J  (4b) 

/ It 

? I'  ‘ 

bj  = [^p  ( p ) - g(  *- )er  f fj  1^1  + 

+ (a, -a  ) ! exp  ( - p ^ ) -er  f c ( i|  )1  (4&) 

1 o i-  ■ J 

where  the  function  g(x)  is  the  solution  of 
the  equation 

p"-t2xp'  - 2erfc(x)  , g(o)=  0 , g’(o)  = l (4  7) 

shown  with  a dashed  line  in  figure  (5). 

We  i.otice  for  future  reference  that 
g(»)  =-  1 . S8. 


The  e .xpr  e s s i on  for  -F^  is  identical  to 
that  of  Fj  with  a^  replaced  by  -a^  . Equa- 
tion (44)  can  now  be  used  to  calculate  the 
unknown  parairioters  a^  , thereby  obtaining 


/v  /(i  + 2A  ) 


= {• 


3/2 


. i ^ (1+  ZA) 

? /7 


-2 


(48) 
(4  9) 


The  asymptotic  expansion  in  Eq.(43) 
can  be  used  to  calculate  large 

times.  The  relation  obtained  wlien  using 
the  first  three  terms  of  the  expansion  has 
been  plotted  in  figure  2,  for  A = 0 , to 
sliow  the  range  of  values  of  e when  it  is 
an  accurate  representation  of  the  exact, 
solut ion . 

A similar  procedure  yii'lds  the  long 
time  mass  fraction  profiles,  written  here 
in  terms  of  the  similarity  variahJe  (i=p/bo 

y - Xga^/irerfc(  M ) + { a^ /iT [^g ( p ) - g ( " )J  + 

/7a  , 

+ /7(aj  + — j--)  erfe  ( p ) ) (l  - “)  + 

^ e 


It  ID  ot:;cr'vc'd  Ut<ii  the  fra-j;ion  { ro- 

f:Te£  dept-nd  orTy  on  tho  jmra;;.oter  A. 

rquolions  (^^0)  and  (^3)  xay  be  usea  to 
derive  an  approxirri.ite  forni  of  the  relation 
valid  for  all  times. 


X'^  + tiY  “b4-d.+d„x 
')?  /XX  \ S 1 J s 

0 = ) 


(51  ) 


w >i  o r e 
a = -J 


2a. 


b = 


(82) 


‘ip  = 


bf  - 2nb  i ag  - aa^ 
o 


lA  p , A.'  ] 

L a a A a ^ - 


(53) 


dj  - d^a  = 


aa  ^ + La  ^ /o 
o o 


7.h/a 


, f-; 


p ^ r Jlil 

a aaZ-' 
o 


( 54) 


f = 


m ^ ( 1 1 4 A / 3 TI  ) 
(2  + A)3 


g - it/(2lA)^  (55) 


Express  ion 
0 t repi'oducc.n 
equat  ioti  (43), 
0 , rejiroduccs 

terms  of  order 


(51)  v;hcn  expanded  foi'  larj*,e 
the  first  t h i' o e terms  of 
•j7id  when  expanded  for  small 
the  correct  expansion  up  to 
X^  . The  approximate  rela- 


tion (51)  when  Civavin  in  fif.urc  (2)  can  net 
be  d i s t i u i s hod  from  the  exact  solution 
obtained  numerically. 


Since  the  obiective  oi’  the  resent 
study  is  the  analysis  of  the  gasification 
process  in  order  to  develop  a follow  up 
in  VC  s t i got  icn  of  gas  ydiase  ignition,  we 
are  specially  interested  in  obtaining  tem- 
perature ai.d  mass  fraction  prof.iics  in  the 
gas  phase. 

Equations  ( 3fc ) and  (41)  valid  for  snail 
values  of  titue  and  cquot.'ons  (4?)  and  (50) 
valid  for  long  times  will  permit  the  analy- 
sis of  eaiTy  ignitions  (very  reactive  pas 
phase)  and  late  ignitions,  thereby  coverinr, 
both  limiting  behoviours. 

Tn  figure  (5)  profiles  of  temperature 
deviations  x » X from  the  inert  or  noii- 
gasification  values  and  the  non-dimensional 
mass  iracticn  are  shown  as  a function 

of  n for  different  values  of  o . We  have 
only  shown  the  case  o = 0 , Lc  - 1 , for 
which  the  equations  and  boundary  corditions 
defining  x i X.,  ? ^nd  y co  i nv'.  i ilo . Using 
figure  (5)  and^oqua t j ons  (22),  (23),  it  is 
easy  to  derive  the  temperature  profiles 
for  a given  3/f  . 

Notice  that  the  fuel  mass  fy*action 
docs  not  t»ocome  of  order  unity  until  very 
long  litres  when  x becomes  oi  order  3 , so 
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ri<*.5  Tompei'-*  ^ I .ir.d  nasr.  fi-actior 
prof  2 2 ' ?'  a - , I.c  = 1 . 


A 0 X p ( • ( i‘  u / r ) ^ ‘ = 

- /■ 


»'  II  u ( 1 - <!  ^ ( T ^ - T ) f c - c ) / L ) / Cl  I’ 
w * o P 


( ■■  ^ 


obtaiiM^U  froM  the  boundary  c ond  i t i o-js  i ri 
E q s . ( 1 V ^ and  ( f 9 ) . 


hqs.(G'’)  and  (02)  cive  A,  B ard  cx 
pllcitely  in  teriiir  of  u,  while  Hq.(09) 
p,ivcs  t(u).  Thus,  we  have  <ar.  explicit 
pc.rr.Tiioiric  repreoontat  ion  of  the  r*'^lation 
invoJvin^  tfiC  parameter  u . Since  u 


is  of  order  unity  when  is  of  order  uni 


ty,  an  explicit  /elation  ^-ivinp  rg(-i),w:t 


eri'ors  of  order  1/3  , cen  bo  obtained  by 
V.' r i t i n g u ~ 1 in  h q . ( 0 9 ) . 


The  fuel  riass  f2’action  is  giver,  by 


* r 


/Y... 


- { ( T)  •-  u ) } 


1 +crf ( u/be ) 


(63 


where  Y.,  is  deterniinod  by 
} s ^ 


Y ~ 1 F s ex p ( - ir  i.e  ) 

u/^r,e  l + erf(u/ho) 


t 6m 


so  that 

(T,-T^),{T^-T^)  (SO) 

is  of  o y,  and  Lhe  solution  given 

by  bqo.(*t.  ■ .ii.  . (SC)  is  lio  lorigor  valid. 

To  obttain  r!»o  soljtion  of  the  system  of 
Eq^.(14)-(21)  for  those  large  values  of 
tijne,  we  write  th?  cysten  in  terms  of  the 
variables  n ~ C/2*^  and  t , and  notice 
iliat  the  time  dei  ivativti  terms  are  of  or- 
der i/s  rel.ative  to  the  remaining  terms. 
Thus,  in  first  approxiiiiat  ion  ve  can  use 
the  quasi-similarity  as-^umptioii  as  wc  did 
v;hcn  founding  the  rsolution  of  the  system 
of  Lqs.(27)-(3  0)  for  large  o. 

The  resulting  system  of  ordinary  dif- 
fciential  equations  can  be  solved  to  yield 


(6,j-G),'r'  - Ae  1 f c ( r.  4 \.i  r.'r)  + erfen 

( b'7  ) 

(0-  -e  ) Berfc(n-u)-erfr. n 

l£  E 

(58) 

Wl.Cl’O 

ur/{ru(l  1 D/t  } = v,'Vj  . 

(59) 

V/Vj  = C- Xp  [ - S /' ^ 1 '1  t 

(60) 

aud  tile  constant  A and  B are  fiwii 
equations 

by  the- 

♦ ^ = Aoi  feC' j/r)4l  [ro.’fc(-u)-l  }/r  (1.1) 


obtained  from  tht-  boundai'V  coiiditicri  in 
Kq.(21).  A siniltU'  equal  .'on  is  obtaio'^d 
to  d e s c r-  i b 0 the  i d i s e r m a s r,  f i'  ^ c t i o ; i b / 

replacing  Y,  and  Y,.  bv  1-Y^/Y,  , and 

1-Y^^,/Yq^  respectively. 

The  time  evolution  of  the  tomnerat arc 
and  mass  fi’ac  lions  at  the  suj'f  4ce  of  a c: 
densed  material  i]j~t  according  f 

an  endothermic  Arrhen'us  dccompos i T ion  n 
action  after  being  exnori.-d  to  a hot  gas, 
has  been  described  analytically  by  r.ieans 
of  asymptotic  techniques. 

Upon  contact  with  the  hot  gas,  the  su 
face  ten.perature  is  raised  to  an  intci-riK.- 
diate  "jump”  value  equal  to  the  one  exist 
ing  in  absence  of  ga s i f £ c a c i on  . Gasifica- 
tion however,  causes  the  surface  tempera- 
ture to  decrease  with  time,  siiico  pa>'t  o 
the  energy  corning  from  the  hot  gas  is  ab- 
sorbed by  tlie  endothern;ic  reaction  and,  i 
addition,  because  the  gas  : f icat  io'n  flow 
bi'ingii  cold  material  towards  the  surfticf 
and  blows  the  hot  gases  away.  lJu>  ing  an 
early  period  in  which  the  decreasing  sur- 
face regression  rate  is  still  o1  the  orb. 
of  its  initial  valuo  at  the  jump  lotrjr.ia- 
turc,  the  Arrlicnius  exponent  can  be  lir.e- 
arixed  around  the  jump  temperature.'  und 

only  the  gradionl  of  the  no  n - ga  s i f i c a t : ci, 
tcmperatuj'e  profile  enters  the  ccaivectivo 
effects.  As  a T'csult  of  these  simpJifire. 
lions,  the  surface  temperature  iiis*ory  ea 
be  described  by  the  solution  of  inte- 
gral equation  riving,  v.  - !(T.-'i'  )RT?  in 


F 


terms  of  a non-d  imetis  iona  1 time  o and  a 

parainoter  A measurini?  the  convective 
effect::;  this  is  re  ^-r  eon  t cd  in  lip,  .2.  Tho 
surfacM;  v,ilue  of  th?  fluid  mass  fracticn 
is  given  in  fig. 3,  Closed  form  expros~ 
sions  of  the  temper-it  ur . and  concentration 
profiles  are  also  given  in  the  form  of  ex- 
pansions foi’  small  and  large  values  of  o . 

When  the  variable  x^.  becomes  large  of 
order  6 > the  complete  convective  ter' ms 
have  to  be  retained  and  the  exact  forni  of 
the  Arrhenius  exponent  has  to  be  used. 
However,  for  those  large  tines  a quasi- 
similar  solution  exist  when  t hj^  distance 
to  the  surface  is  scaled  by  ?/i  , so  tliat 

explicit  expressions  have  been  obtained 
for  the  t ‘j rr.[’e r a t ure  and  coiic e n t ra  t i o n pro- 
files in  tei'ins  of  the  similarity  variable 
n and  a function  u of  time  given  by  Eq. 
(59). 

Ve  consider  in  an  Appendix  the  case 
when  tho  gasificatic)n  5-^  assumed  to  obey 
an  equilibrium  C 1 a s i u s - C 1 a pe  yi'on  condition 
instead  of  the  regression  rate  law  of  Eq. 
(9).  The*  t .•'mperature  and  mass  fraction 
profiles  arc  described  in  terms  of  a sim- 
ilarity solution  valid  for  all  times,  _ 

involving  the  similarity  variable  n = ^/2/t. 

Tlie  surface  temperature  does  not  change 
with  time  after  jumping  at  t = 0 to  a val- 
ue determined  by  gasification  effects.  The 
regression  rate  is  found  to  vary  as  1//t"  . 
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A p 0 n d 

In  this  appendix  vk-  present  an  anal- 
ysis  of  the  cases  whu.  phase  equilibrium 
is  ina  i a t a i n 0 d du!‘ing  the  gasification  pro£ 
ess.  In  practice  the  equilibrium  analysis 
will  be  usually  applicable  to  the  case  of 
li^^uid  fuels  end  the  rate  process  analysis 
will  he  .better  suited  f o i’  the  c a 5; o of 
solid  propellants,  plastics  and  other  sol- 
id mat  cr  ial s . 


_P, o , 

T - 'I 
go  o 


1 - C ^ e r'  f c ( n - V ) 


e r f c ri  - V ) 

y = 

crfc(-V /l? ) 


and  the  boundary  conditions,  yield 


bs  = r ^ 


exp( -V^  Le  ) 

v/tTlc  erfc('V/  L^e  ) 


T - T 
2 o 


C,crfc(— ] = l-C^erfc(-V)  = 

j ^ r 2 


Pe  /erfeC-V)-  VFt'V/ra' 
fe  ^ /erfc(-V)fexp{-(rv/r)^j/erfc(rv/r) 


The  C.las  i us-Clapcyrcn  equation  relates 
the  svrf.ace  mass  ti*a''tion  to  the  surface 
tcmpvr nture  through 


c(T  -T  ) 
go  o 


where  wc  have  made  the  non-essential  assump 
tion  that  the  molecular  weights  of  the  ga£ 
ecus  species  are  tqual.  If  in  addition  L 
changes  slowly  in  the  range  of  integration, 
then 


= exp[L/P,Tj,-  I./rTj 


Equation  remain  va.lid,  but  equa- 

tion (6)  is  I'cplaced  by 


V \ £2.  I 

S E 3>l'  I , 


LC  3Y_| 


This  problem  is  simpler  than  the  one 
in  the  main  text  since  it  admits  a similar 
ity  solution  valid  for  all  times.  It  turns 
out  that  the  rcgi’osricn  velocity  varies 
with  time  as  \ f and  the  temperature  and 
concenti’at  ion  profiles  become  I : mc  inde- 
pendent if  described  in  terms  of  the  simi- 
larity variable  p = 5/2/r,  Introducing 


From  equation  (A9),  the  value  of  1’^/1’q 
is  obtained  as  a function  of  V for  fixed 
values  of  ci',r,F  • Equation  (A2)  is  then 
used  to  deduct’  To/T  as  a function  of  e’, 

D O 

r , r , Le  , L/RT  , and  T /I  . Once  this 


r , r , Le  , L/RT  , and  T /I  . Once  this 
* ’ * o go  o 

curve  is  known,  the  value  of  V for  the 

given  Tp  is  calculated,  and  therefore  the 

value  of  'T'  /T  can  be  deduced  from  equa- 
s o 

ticn  (A9).  By  this  procedure  the  tempera- 
ture and  iiiass  fraction  d i s t r ii.'u  t ion  aic 
knov/n  fer  all  time. 

Since  the  solution  to  the  g,a  s i f i c a t ion 
problem  involves  several  parameters,  ther£ 
by  possessing  many  interesting  limiting 
behaviours  whose  analysis  will  considera- 
bly lengthen  the  present  paper,  a para- 
metric study  of  gasification  under  equilib 
rium  conditions  is  postponed  to  a subse- 
quent publication. 


3Y,/3nl. 


the  solution  to  equt.tions  (lV(d)  is 
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